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Abstract: We consider a three-dimensional effective theory of Polyakov lines derived 
previously from lattice Yang-Mills theory and QCD by means of a resummed strong coup¬ 
ling expansion. The effective theory is useful for investigations of the phase structure, 
with a sign problem mild enough to allow simulations also at finite density. In this work 
we present a numerical method to determine improved values for the effective couplings 
directly from correlators of 4d Yang-Mills theory. For values of the gauge coupling up to 
the vicinity of the phase transition, the dominant short range effective coupling are well 
described by their corresponding strong coupling series. We provide numerical results also 
for the longer range interactions, Polyakov lines in higher representations as well as four- 
point interactions, and discuss the growing significance of non-local contributions as the 
lattice gets finer. Within this approach the critical Yang-Mills coupling f} c is reproduced 
to better than one percent from a one-coupling effective theory on N r = 4 lattices while 
up to five couplings are needed on N r = 8 for the same accuracy. 
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1 Introduction 

Yang-Mills theory as well as QCD can be mapped to three-dimensional effective Polyakov 
loop theories by integrating over the spatial gluon degrees of freedom, both in the continuum 
[1-4] and on the lattice [5-10]. The resulting effective theories describe the physics of the 
phase transitions and thermodynamic behaviour. They are simpler in the sense that their 
non-perturbative treatment requires less numerical effort and, at least in some limiting 
cases, the effective interactions can be associated with the low energy degrees of freedom. 
Since part of the degrees of freedom have been integrated out the sign problem, which 
hampers lattice simulations at finite baryon density, is much milder in the effective theory. 
Indeed, at least for very heavy quarks [8] and in the chiral limit on very coarse lattices 
[11] analytically derived effective lattice theories allow for simulations of the nuclear liquid 
gas transition in the cold and dense regime of QCD. On the other hand, effective actions 
generically include infinitely many and arbitrarily involved interaction terms, restricted 
only by the symmetries of the underlying theory, and become practical only after some 
drastic truncation. It is therefore essential to understand the relevance of the various 
terms and the effects of truncations. 

In this work we consider Yang-Mills theory on a lattice with periodic boundary condi¬ 
tions and the effective action obtained by integrating over the spatial links, 

e -S eS [W] = J [dU . ]e -S[u\ ; (Li) 
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where S[U] is the standard Wilson lattice gauge action and depends only on the 
Polyakov lines W(x) = n yu - y=x Uo(y). It contains all n-point interactions at all distances 
and with all powers of loops that are consistent with the Z(iV)-centre symmetry of the 
Yang-Mills theory. The expectation value of all observables represented by Polyakov lines 
can be calculated in the effective theory, 

(0) = ^ J[dW]e~ s ^ w ^O{W) = ^ J[dUo}[dUi]e- s ^O(W[U]) . (1.2) 

The normalisation of the path integral is chosen such that 

Z = J[dUi\e- s M = J[dW]e~ s ^ w ^ . (1.3) 

The constant part of the effective action, which does not depend on the W. can be neglected 
(absorbed in the normalisation) if one is interested in the phase transition or VP-dependent 
observables only. However, for detailed thermodynamic properties the T dependence of 
the constant part might be relevant. 

Svetitsky and Yaffe conjectured that the short range Polyakov interactions are the rel¬ 
evant terms for the phase transition [5]. This is based on the fact that these are the 
dominant contributions both in the strong and weak coupling limit, while for all couplings 
interactions are screened by the mass gap of the theory (see also the discussion in [12]). 
In [7] the integration of the spatial links was performed by means of a strong coupling 
expansion, where interaction terms are parametrically suppressed with increasing distance 
and n-points. Simulation of the effective action truncated to the leading nearest neigh¬ 
bour coupling gives already a good description of thermodynamic functions and the phase 
transition, allowing for a continuum extrapolation of the critical temperature. This already 
vastly improves over direct strong coupling calculations of thermodynamic observables [13]. 
On the other hand, the missing higher order couplings and long range interactions hamper 
the applicability of the effective theory to correlation functions and the string tension, 
where non-local contributions play an increasing role as lattices get finer [12]. The same 
observation was made with a non-perturbatively determined form of the effective action 
[9, 10]. 

In this work we present a new method to numerically determine the couplings of the 
effective theory directly from simulations of the 4d Yang-Mills theory. The method employs 
the fact that there are different character expansions of the effective action. In one of them 
the effective couplings emerge as expectation values of n-point functions of Polyakov loops 
in the full theory, which can be measured easily. The relation between effective couplings 
corresponding to different character expansion schemes is then established perturbatively in 
terms of a rapidly converging power series. This allows us to check the range of applicability 
of the strong coupling approach [7, 12] and significantly improve on the results. 

This paper is organised as follows. In section 2 we discuss the character expansions 
pertaining to the strong coupling approach and the new numerical method. We invest¬ 
igate two-point interactions of Polyakov loops in the fundamental and anti-fundamental 
representation in section 3 . The new method is introduced for the one-coupling effective 
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theory and then generalised to include two-point couplings over larger distances. In Sec¬ 
tion 4 we discuss the significance of higher representations as well as three- and four-point 
interactions in the fundamental representation before concluding in section 5. 

2 The different character expansions of the effective action 

The effective action can be expanded in terms of the characters \r of all irreducible rep¬ 
resentations of the gauge group. The irreducible representations of SU(N) are labelled by 
a vector r of N-l canonical labels. In SU(3) one has two-dimensional vectors like (1, 0) and 
(0,1) for fundamental and anti-fundamental representation or (1,1) and (2,0) for the ad¬ 
joint and sextet representation. The characters of the fundamental and anti-fundamental 
representation are the ordinary Polyakov loop T(x) = X(i,o)(f^( x )) = tr[W(x)] and its 
complex conjugate U = X^foy(W) = X(o,i)(W)- Since products of the characters of the 
irreducible representations can be transformed into sums, the effective action can be trans¬ 
formed into a representation without products of \ r (IT(x)) at the same lattice point x. 
The effective action generically includes arbitrary interactions of characters at all distances. 
A possible representation of the effective action that includes all of these terms is thus 

n 

S eS = Y,Y, J2 C [n,xi],...,[r, I ,x n ]Xr(bP(x)) nXr ; (hP(x + Xi)) . (2.1) 

x,r n i= 1 

[ r l ,x i],[ r n,Xn] 

The coefficients of this expansion are not completely independent but related by the sym¬ 
metries of the theory. In our case, the Z(N) centre symmetry restricts the independent 
combinations of the r l; and the cubic symmetry (the remnant of rotational and translational 
invariance in the continuum) restricts the combinations of x,. 

Our aim is a numerical determination of the coefficients in the effective action from 
expectation values in the full theory. The relations between the measured expectation 
values and the effective action should be as simple as possible. Certain observables like 
the Polyakov loop correlators can be obtained with a high precision. The form (2.1) is, 
however, not well suited for such an approach. 

Alternatively, also the exponential of the effective action can be represented by inter¬ 
action terms of characters. This corresponds to a Taylor expansion of the right hand side 
of (1.1), where products of characters at the same point are converted into sums of single 
characters. In this way the expansion of the effective action (2.1) is converted into an 
expansion of its exponential, 

e- S - = A„(l + j:5: ]T A[ xi , ri] . K , rn] Xr(ir(x))nXr.(ir(x + X.;))) . (2.2) 

' x,r n i= 1 ' 

[ r i. x i]> 

■ ■ • j [ r n ; x n] 

The sum is performed in such a way that each combination of characters and lattice 
points appears only once. Because of the orthogonality of the characters, expectation values 
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of the characters in the full theory now directly project on the corresponding coupling 
constant, 


(XrWOjjllXr^Xi))) 

2=1 


1 

z 

1 

z 


J[dU^e-^xAWm fl Xr t m^)) 


(2.3) 


i=l 

n 


e- S ^Xr(W( 0)) II Xri(W(Xi)) = A[ Xi)ri]r .. :[Xn)rn] 


i =1 


(The constant part Ao = /[dt/^]e — has been factored out to arrive at the proper norm¬ 
alisation of the expectation values and will not be considered in the further analysis). The 
disadvantage of this representation of the effective action is, however, that any truncated 
version of the expansion is not appropriate to reproduce the phase transition. The range 
of interactions in the corresponding effective theory is restricted to short ranges. At larger 
distances all correlations, in particular the Polyakov loop correlator (L(O)iA(x)), are zero. 

A different representation, which includes all long range correlations already in the 
truncation to its first term, follows in a natural way from the strong coupling approach [7] , 


o-Seff _ 


x,r,n 


= AoII II ( 1 + A [ Xl ,r 1 ],...,[x n ,r n ] Xr(VP(x))[]x ri (lT(x + X i ))+C.C 


[ F 1 ’ X ll ’ 


2=1 


...,[r n ,x n ] 

(2.4) 


The product J} now runs over all inequivalent combinations of irreducible representations 
r and lattice positions x. The coefficient of the trivial representation has been factored 
out since it is irrelevant in the calculation of expectation values. A truncation in this 
representation is a restriction of n and r,ri,...,r n to a certain number of irreducible 
representations and of the x,; to a finite number of interaction distances. Such truncations 
arise naturally in a strong coupling expansion. This representation combines the advantages 
of (2.1) and (2.3): long range correlations are included and the calculation of the coefficients 
from expectation values is still feasible. The Polyakov loop correlator in such an effective 
theory is non-zero at arbitrary distances and decays exponentially, as in the full gauge 
theory. 

All discussed representations are complete expansions of the effective action, implying 
relations between the coefficients c, A, and A. If the expansions are truncated, only approx¬ 
imate relations can be established. In this work we use the representation (2.3) in order to 
determine a finite number of the couplings A numerically by simulating the corresponding 
n-point functions. In a second step, we evaluate the same expectation values as a perturb¬ 
ative series in small A using the representation (2.4). The relations A(A) are then realised 
as power series which can be easily solved for the A. In this approach the effective coup¬ 
lings, containing information of the full correlators, clearly deviate from those obtained in 
the strong coupling expansion. But we find them to be sufficiently small over the entire 
parameter range of /3, N T to justify a perturbative expansion, with an ordering essentially 
the same as in the strong coupling expansion. In particular, the dominant contribution is 
the nearest neighbour interaction. 
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Above the phase transition strong coupling series for the effective couplings are no 
longer valid, whereas their numerical determination is of course still possible. However, 
the effective action now has to take into account the non-zero expectation value of the 
fundamental loop (as well as that of Polyakov loops in higher representations with non¬ 
zero N-ality). One possibility is to add a one-point term, here we include the expectation 
value by a simple shift of the Polyakov loop variables, 

L -> L - (L) . (2.5) 

It is important to note that the effective theory (2.4) captures the signal of the phase 
transition even without this shift, the non-vanishing expectation value resulting from the 
non-perturbative dynamics. By contrast, this is impossible in the representation (2.2) 
without addition of a one-point term. Our main interest is in the effective theory up to 
the phase transition and its prediction for the transition point. For a better presentation 
of the effective couplings we include in some cases data above the transition, including the 
shift (2.5) . 


3 Two point interactions in the fundamental representation 

The simplest contributions to the effective action are two-point interaction terms of Polyakov 
loops in the fundamental and anti-fundamental representation at all distances. We first 
list the strong coupling results and the corresponding ordering of the couplings, then we 
compare it with the non-perturbative numerical data. 


3.1 Strong coupling results 

The expansion parameter of the strong coupling series is u(/3) = ^ + ... < 1, the coefficient 
of the fundamental representation in a character expansion of e ~ Sp ^ with the plaquette 
action S P (U ) [14]. To leading order the nearest neighbour coupling is u N t , where N T is the 
number of lattice points in temporal direction. The leading order contribution to couplings 
over larger distances is at least dN T , where d is the shortest distance along the links of the 
lattice (“taxi driver distance”). We write the effective couplings always as a product of 
leading order contribution and corrections. For the nearest neighbour interaction we have 
computed the corrections up to high orders in u, they can be resummed into an exponential 
factor (for details see [7]), 


X = xT°) 

Al - '"'[(1,0,0),(0,1)] 


4 

u exp 


N T 

u T exp 


4(4 u 4 + 12m 5 - 14m 6 
W(4m 4 + 12m 5 - 14m 6 


7 a 

36m 7 + —m 8 + 

7 295 o 

36m 7 + —m 8 + 


1851 

~To 

1851 

U 


-XT + 


-M y + 


1035317 

5120 

1055797 

5120 



N t = 4 

N t > 6 . 
(3.1) 


Here we introduced a simplified one-index notation for the couplings A, which we also apply 
for the representation (2.2) of the effective action (Ai = Ar^° 0 0 j ^ etc.). 
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The next terms are interactions at distance y/2 and 2, which have the same “taxi-driver 


distance” d = 2: 

f tt 10 (12 + 26u 2 + 364u 4 ) : N r = 4 

A 2 = A j^o),(o,i)) = u u (30 + 66n 2 ) : N r = 6 (3.2) 

[ N t (N t - 1 )u 2Nt+2 : N T > 6 

As = A[ l£l m0tl)] = 4AU 27V ” +6 (1 + 12 u 2 + (8 JV t + y)u 4 ) . (3.3) 

We also consider interactions with d = 3 

A 4 = A = u 3N ^ +4 {N t 4 + 6 N T 3 - 13 N t 2 + 6N t ) (3.4) 

As = A[(2^1,0),(0,1)] = N r( N r ~ 1)^ +4 (3-5) 

Ae = a[J’ 2 >0)>(0i1)] = 4N t u 3Nt+8 (1 + 12u 2 ) . (3.6) 


Note that the corrections to the leading order strong coupling result have been computed 
with different precision for the various couplings since their calculation gets much more 
involved at larger distances. In the following we limit our investigations to a maximum 
order of 41V T on the N r = 4 lattice, which excludes the interaction Ag. The effective action 
(2.4) has thus the following form 

e _ ‘ S eff ) = ^1 + Ai (L(x)L(x + i)^-h L(x)^L(x + i))J 

x,i= 1,...3 

II (l + ML(x)L(y? + L(x)iL(y))) 

[x,y] 

(l + A 3 (L(x)L(x + 2i)t + L(x)^L(x + 2*))) 

x,i= 1,...3 

n (1 + A 4 {L{x)L[y)^ + L{x) ] L{y))^j 

{x,y} 

II (l + HL(x)L(y)' + L(x)'L(y)j) , (3.7) 

<x,y> 

where [x, y\ indicates the product of all points with distance \/2 (i. e. connected around 
one corner), the product {x,y} contains all points with distance \/3, and < x, y > with 
distance y/5. 

3.2 The effective coupling in the one-coupling action 

The simplest truncation of the effective action contains only the nearest neighbour coupling 
Ai, setting all other couplings in (3.7) to zero. Previous investigations based on the strong 
coupling series (3.1) have shown that this simple form of the action allows already a de¬ 
scription of the phase transition to an accuracy of about 10% for lattices with N T = 2 — 12 

[7]- 

On the other hand, the exact numerical value of the coupling in the expansion (2.2) is 
determined by the simple correlator, 

A 1 = (L((0,0,0))L((l,0,0))t), (3.8) 
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Figure 1. Numerical solution for Ai in the one-coupling effective theory compared to its strong 
coupling series (3.1) for N T = 4. Different truncations of (3.9) and inclusion of some larger distance 
couplings (3.11) (including X n ) are shown. The horizontal line marks the critical effective coupling, 
vertical lines the critical /3 C in full Yang-Mills theory and (3^ predicted by A^ [7]. 

which we measure in a simulation of the full 4d theory. The same correlator can be 
calculated in the effective theory up to O(Aj) in an expansion of X±. Equating these 
two expressions one arrives at the following relation between the coupling constants in 
representation (2.2) and (2.4), 

Ai = Ai + 4Af + 24Af + O(A^) , (3-9) 

which can be solved numerically for X±. Note that Ai is then exact (up to statistical 
errors and finite volume corrections), whereas a systematic error is introduced to Ai by the 
truncation of the right hand side. 

Figure 1 shows the results obtained for an N r = 4 lattice. The strong coupling series 
gives a reasonable approximation to Ai in a wide parameter range. We observe that our 
numerical approach is self-consistent: even though Ai deviates from the strong coupling 
expression (3.1) at larger u, its value remains small and the truncation of (3.9) is justified. 
Orders beyond O(Af) in (3.9) are not relevant and the influence of the interactions at 
larger distances (A n ) appears to be suppressed. Details of the computation including these 
interactions are explained in the next section. 

The phase transition in the one-coupling effective theory is at (Ai) c = 0.1879 [7], i.e. the 
intersection of the curve Ai(/3) with this value gives the prediction of the one-coupling 
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theory for f3 c . Note that for the numerically determined coupling Ai this coincides with 
the value at the true f3 c of full Yang-Mills, Ai(/3 C ) ~ (Ai) c . Thus, the one-coupling model 
with the effective coupling determined by our numerical approach provides a quantitatively 
correct prediction for the transition point at N T = 4 and significantly improves over the 
strong coupling result, which tends to overestimate j3 c . However, we shall now see that for 
larger N T , i.e. on finer lattices, more couplings are necessary for the same precision. 

3.3 Interaction terms at larger distances 

We now include all two-point interactions of fundamental and anti-fundamental loops whose 
strong coupling series contain contributions up to 0(u 4Nt ), (3.7). Again we directly determ¬ 
ine the couplings A n in representation (2.2) by the corresponding Polyakov loop correlators, 

A 2 = (L(( 0 , 0 , 0 ))L((l,l, 0 ))t); A 3 = (T(( 0 , 0 , 0 ))L(( 2 , 0 , 0 ))t) .... (3.10) 

As in the one-coupling case, the same correlators are computed in the effective theory 
(3.7) using an expansion in the effective coupling constants A n . However, in this case the 
truncation is less obvious. With infinitely many couplings there are also infinitely many 
mixed terms with the same overall number of A’s, i.e. there is no simple power counting. 
We therefore use again the strong coupling expansion as a guiding principle and truncate 
terms whose leading order contribution would exceed 0(u 4Nt ). In that way a coupled 
system of equations is obtained, 

Ai = Ai + 4Af + 8 AiA 2 + ... 

A 2 = A 2 + 2\j + 4A| + I 6 A 4 + 36A?A 2 + ... 

A 3 = A 3 + Af + 4A| + 24AfA 2 + 12\f + ... (3.11) 

A 4 = A 4 -I- 6 Af T 6 AiA 2 T ... 

A 5 = A 5 + 3Af + 2 AiA 2 + ... , 

that can be solved for the couplings A n . To estimate the effects of the truncation in the 
expansion of the couplings we have also considered a second set of equation, 

Ai = Ai + 4 A 4 + 8 AiA 2 + ... 

A 2 — 2 A 4 = A 2 + 4A^ + 4A^A 2 + ... 

A 3 — A 4 = A 3 + 4A^ + 8 A^A 2 + 4 A 4 + ... (3.12) 

A 4 — 6 A 4 = A 4 + 6 AiA 2 + ... 

A 5 — 3Af = A 5 + 2AiA 2 + ... , 

where the complete nearest neighbour correlation is subtracted. This corresponds to a 
resummation of certain dominant higher order terms. We get consistent results for both 
sets of equations, which is an indication that the truncation errors of the system of equations 
is small enough. In the following we will use (3.11) for the determination of the A n . Figure 
(2(a)) shows the corresponding results and illustrates the accuracy of the Polyakov loop 
correlators required in order to compare with the the strong coupling predictions. The 
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(a) Measured correlations (N T = 4) 


(b) Effective couplings (JV T = 4) 



(c) Logarithmic representation ( N r = 4) (d) Effective couplings (N T = 6) 


Figure 2. (a) Effective couplings defined by correlators, (2.3), from simulations of 4d Yang-Mills 
on N t = 4. (b) Solution of equations (3.11). Lines show the strong coupling series for the couplings, 
(c) Same as (b) in logarithmic representation, (d) Same as (b) for N r = 6. Vertical lines mark the 
critical coupling u(/3 c ) in full Yang-Mills theory. 

numerical solution for the set of equations (3.11) is shown in Figure 2(b) and 2(c). At small 
j3 the result agrees with the prediction from the strong coupling expansion, which provides 
a good approximation for the couplings Ai and A 2 over a wide range. The disagreement 
between the numerical determination and the strong coupling expansion is larger for the 
long range interactions and presumably due to the much shorter series for these couplings. 

While the phase transition on an N r = 4 lattice is predicted with a good accuracy in the 
one-coupling effective theory, the long range interactions are necessary to reproduce the 
Polyakov line correlator of the 4d Yang-Mills theory. This is illustrated in Figure 3. In [12] 
we observed that the mis-match of the correlator remains when additional couplings are 
added in the strong coupling approach, since their leading order contributions are too small. 
The situation is different with the numerically determined effective couplings. While the 
correlator in the one-coupling theory still shows large deviations from the full one, inclusion 


- 9 - 


























1 


0.1 “ 



0.01 

g 

0.001 

g 


o' 

0.0001 

g 

le-05 


le-06 


le-07 


0 


Ai strong coupling approach 

A] 

Al,A2 

Al,A2,A3,A4 

full YM 


+ 

* 


f ! 


t 




t . 


¥4 x 


I f ° 


- 

-G ■ ■ 


9 

1 □ 





-i. 

+ 

■ i 9 < 

>1 

T 1 



* i : 

Mid 


- 1 -' 


i 

5 

* 

ft 


3jf- 

| < 

- 

) i 




* ! 

: i 


9 



R 


Figure 3. Polyakov loop correlator in full Yang-Mills and the effective theory (3.7) at (3 = 5.4 
on a 24 3 x 4 lattice. The different truncations of the effective theories are based on the couplings 
determined by (3.11), as well as the strong coupling series (3.9) for the one-coupling theory. 

of two-point interactions up to A 4 results in good agreement. Thus the main improvement 
achieved by our numerical procedure are better estimates of the longer range interactions, 
where no sufficiently long strong coupling series are available. 

3.4 Different N T and the continuum limit 

Before we turn to the more involved interaction terms, we discuss the behaviour of the 
effective couplings towards the continuum limit of the 4d Yang-Mills theory. Since the 
strong coupling series for larger distance interactions are parametrically suppressed by 
additional powers of N r compared to the nearest neighbour interaction, it is tempting to 
expect their relevance to diminish with growing N T [7]. However, Figure 2 shows that at 
the phase transition all couplings become of comparable size. This is not unexpected, since 
the critical coupling represents the convergence radius of the strong coupling expansion. 
The general behaviour for growing N r can be inferred by inspecting (3.11). The A* on the 
left side of the equations are given by bare Polyakov loop correlators. For fixed temperature 
in the 4d Yang-Mills theory, these get smaller with growing N T , i.e. the continuum limit of 
(3.11) is zero. The numerical values of the couplings for N T = 6 are shown in Figure 2(d). 
The effective couplings at a fixed u are smaller than for N T = 4, a behaviour also shared 
by the corresponding strong coupling series. Consequently, all couplings remain small on 
the way to the continuum. 
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Figure 4. Numerically determined (Ai(u(/?)), A 2 (it(/3))) for various N T . The grey line repres¬ 
ents the phase boundary of the two-coupling effective theory. Red points mark the critical region 
0.995 u(f3 c ) < u < 1.005m(/ 3 c ) in full Yang-Mills theory for the different N T . 

Consider now the one-coupling theory for growing N r . Because Ai is diminishing with 
N t for fixed physical T, equation (3.9) will at some point no longer have a solution corres¬ 
ponding to the critical value (Ai) c for the phase transition in the one-coupling theory, hence 
a second coupling has to be added. The two-coupling effective theory has a critical line as 
shown in Figure 4. The intersection of the curves Ai(A 2 ) for given N r now give a prediction 
for /3 C in implicit form, which allows (Ai) c to shrink with N T , but at the same time (A 2 ) c 
has to increase. This is remedied by adding a third coupling, etc. The general situation is 
thus: any effective theory with a finite number of couplings features a critical hypersurface 
representing the phase boundary between an ordered and a disordered phase. For fixed 
N t , the effective couplings A i intersect these hypersurfaces to provide the predictions for 
the critical couplings. Towards the continuum limit all intersection points have to move 
towards the origin, which enforces the addition of more couplings as N T is increased. 

How many couplings are required in practice depends on the demanded accuracy. In 
Figure 5(a) the critical coupling for N r = 6 with different truncations of the effective 
action is indicated by the peak of the susceptibility. The transition point with two and 
more effective couplings is within 1% of the full theory, while the truncation to one coupling 
does not show any transition in the considered range. The importance of interaction terms 
at larger distances further increases at N r = 8, Figure 5(b). Agreement to 1% between 
effective and full theory is obtained by including the couplings up to A 4 . This is also 
consistent with Figure 4. The effective couplings obtained in the 5% critical region of 
Yang-Mills theory nicely coincide with the transition line in the two coupling model for 
N r = 4 and N T = 6. For N T = 8 the truncation to the two-coupling effective action is not 
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Figure 5. Susceptibility of the Polyakov loop. The x-axis gives the relative deviation from (/3 C ) of 
the full Yang-Mills theory. 




Figure 6. Red points give the numerical solution and black points the strong coupling series for 
Ai (/3 C (N T )) . The shaded areas give the region of (3 within 8% (left) or 5% (right) of the true critical 
value /3 C . 

enough for an accurate description of the phase transition. 

Finally we use the numerical coupling Ai from our least truncated effective theory to 
assess the quality of the one-coupling theory with the analytic expression for Ai. Figure 
6 compares the numerical and analytical mappings, Ai (j3 c (N T )), both evaluated at the 
true /3 c (N r ) of the full theory. Interestingly, the quality of the analytic function appears 
to initially improve with N r , until it breaks down for N r > 10, where Ai(/3 C ) ceases to 
get smaller. Also shown are regions Ap(/3 C ± A), where A amounts to 8% (left) and 5% 
(right) deviation from the true value. The horizontal line is the critical value (Ai) c of the 
one-coupling theory. Hence, permitting 8% relative error in the prediction for f3 c , the one- 
coupling theory with the analytic expression for Ai may be used 1 , whereas more couplings 
and longer series are required for more stringent accuracy requirements. 


1 Indeed, the predictions for /3 C from the one-coupling theory in [7] are all within 7-8% of the true result. 
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Figure 7. Left: Correction P c (u), (3.16), to the strong coupling result P sc (u) of the parametrisation 
(3.13) in a logarithmic representation. For each N T the curves are restricted to u < 0.995it(/3 c ). 
The curves contain the strong coupling results to order u 10 and fits of the orders u 11 to u 24 . Right: 
Deviation from the fitted polynomials in the critical region. The simulations were done on 4 x 24 3 , 
6 x 24 3 , 8 x 24 3 , 10 x 30 3 , 12 x 36 3 , and, to check the finite volume effects, 6 x 36 3 , 8 x 48 3 . 

3.5 Functional form of the nearest neighbour coupling 

After resummation of the strong coupling series, the nearest neighbour coupling in equation 
(3.1) is parametrised as 


Ai(ti, N r ) = u Nt exp ( N r P(u , N r )) , (3.13) 

with P(u,N t ) a polynomial in u for every value of N r . Between N T = 4 and N T = 6 
the coefficients of these polynomials differ only at order u 10 . For larger N T the difference 
is at even higher order such that P(u, N T ) ~ P(u) can be considered as approximately 
independent of N r , as long as the expansion and resummation converge. 

In order to check the validity of the parametrisation (3.13), consider the coupling Ai 
corresponding to the Polyakov line correlator, which is related to the renormalised free 
energy of a static quark anti-quark pair at a distance R of one lattice spacing a, 

Ai = (L((0,0,0))L((l,0,0))t) = Z(u) n t exp(-F R (R = a,T)/T) , (3.14) 

with Z(u) the multiplicative renormalisation factor. To leading order Ai = Ai, which 
implies 

P(u) « In (Z(u)/u) — o,Fr(R = a,T) . (3.15) 

The parametrisation (3.13) with a lV r -independent polynomial is hence appropriate as long 
as ciFr(R = o, T) does not show a strong N T dependence, which is generally true for the 
Coulomb part of the short distance region, RT <C 1. 

To investigate the validity of the parametrisation quantitatively, we calculate the cor¬ 
rection P c to the polynomial P sc of the strong coupling approach, 

P(u) = P sc (u ) + P c (u) , (3.16) 
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from the numerically determined P(u). P c (u) is shown in Figure 7 (left) based on sim¬ 
ulations with N t = 4, 6, 8, 10, and 12 for u < 0.995 u(/3 c (N T )). In this entire range the 
data for P(u ) are indeed IV T -independent to a very good approximation. In the immedi¬ 
ate vicinity of the phase transition the parametrisation (3.13) breaks down, as the zoom 
into this region in Figure 7 (right) shows. This IV T -dependence is associated with the T 
dependence of clFr(R = a, T), which in the transition region gets more pronounced due 
to the increasing influence of the deconfined phase. The strong coupling series has only 
information from the confined phase and hence underestimates this effect near the radius 
of convergence. 

Figure 7 (right) illustrates that the -/V T -dependence of P(u) is reduced with increasing 
N t . This also follows from (3.14): the larger N r , the smaller is the distance R = a 
in physical units and the short range part of the free energy is dominated by the T- 
independent Coulomb contributions. At the same time the interactions at larger distances 
become increasingly relevant at the transition point. 

The parametrisation by P{u) serves to extract the predictions for the critical coupling 
(3 C from the one-coupling effective theory as in [7]. As long as its highest order coefficient 
is positive, there is always a /3 C for which 

n ^J\i{u{P c ),N t ) = uexp(P(u)) = N \j (Ai) c . (3.17) 

Conversely, (3.17) with the fitted polynomial P(u) can also be used to approximately 
compute X\(u(P), N t ), assuming that the one-coupling effective action is enough for a 
reasonable approximation of the full theory. 

Using the parametrisation (3.13), the -/V r -dependence can be factored out and the predic¬ 
tions for Ai at different N r can be combined as in Figure 8. Close to the phase transition 
the numerically evaluated effective coupling is smaller than the values suggested by the 
matching of the phase transitions of the full and the effective one-coupling theory using 
(3.17). This deviation is due to the neglected -ZV T -dependence of the parametrisation (3.13) 
close to the phase transition, cf. Figure 7 (right), as well as the increasing relevance of the 
interactions at larger distances. The analytic strong coupling prediction shows the correct 
dependence of Ai on u and N r for small u. At larger u it first underestimates and then 
overestimates the slope of the numerical curve. Due to this fact there is an intersection 
point between the expected transition line from the matching of the full and effective the¬ 
ory and, around N r = 14, the one-coupling action with effective coupling of the strong 
coupling approach reproduces the critical coupling of the full Yang-Mills theory. In the 
region of N r = 8 to N T = 14 the full and the effective theory data seem to converge. 

By contrast, the numerically determined effective coupling suggests that higher order 
coefficients of P(u) are negative, as for P sc truncated at order u 7 . Hence a simple continu¬ 
ation of P(u) does not provide an intersection point with condition (3.17) and the effective 
theory truncated to the nearest neighbour interaction fails to reproduce the transition of 
the underlying full theory at larger N r , as discussed in the previous section. 
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Figure 8. Same data for Ai as Figure 7. In addition, the effective coupling determined by a 
matching of the transitions in the effective and full theory is shown (YM f3 c ). These values are a 
solution of (3.17) for Ai with /3 C determined in simulations of the full theory [15-17]. The leading N T 
dependence has been removed assuming (3.13). The horizontal line corresponds to the transition 
in the effective one-coupling theory for N T = 14. 

4 Higher representations and n-point interactions 

An effective action in terms of two-point interactions seems to provide reasonable results for 
the discussed correlation functions and the phase transition for N r = 4 and 6. On the other 
hand, when interactions over two and more lattice spacings are included, the restriction to 
two-point functions is no longer justified within in the strong coupling expansion. Terms 
involving higher n-point couplings appear at the same order as the long range interactions. 
In the following we investigate the contributions of higher representations as well as n-point 
interaction terms. 

4.1 Strong coupling results 

The n-point interaction with n > 2 and higher representations of the Polyakov loop are 
parametrically suppressed compared to the two-point nearest neighbour fundamental inter¬ 
action. However, relevant contributions appear already at the same order as the coupling 
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A 2 , such as the two-point interactions of higher representations, 


\i,i) = ^ 


- x (i,i) 


[( 1 , 0 , 0 ),( 1 , 1 )] 


= v n ^1 + 9 -N t -), 


^( 2 , 0 ) = ^ 


- x(2,0) 


[( 1 , 0 , 0 ),( 0 , 2 )] 


= w Nt 1 + 6 N t — 
w 


where v = ^(w 2 — u 3 ) + ^-rt 4 + 0(u 5 ) 

8 32 

. 3 2 9 4 297 5 6 

and w = -u 2 - —u 4 + —u 5 + 0{u 6 ) 
4 16 80 


(4.1) 

(4.2) 

(4.3) 

(4.4) 


Polyakov loops of higher representations with n-ality zero can have non-vanishing expect¬ 
ation values in the confined phase. These terms appear in the effective theory as one-point 
couplings. The simplest example is the one-point term of the adjoint loop, 


A (u) = 108 N t {N t + 1 )u 4Nt+2 . 


(4.5) 


Since it is of high order we neglect it in our considerations. 

In addition, there are higher n-point interactions in the fundamental representation. 
The leading three point interaction is a combination of fundamental, anti-fundamental, 
and adjoint loop: 


A (1,0) - —AN „21V T +8 

A [(1,0,0),(1,!)],[(!, 1,0),(0,1)] _ /±JV T« 


(4.6) 


This is beyond the order considered in our current investigations. Four-point interactions 
can appear in terms of two fundamental and two anti-fundamental loops. In the leading 
contribution these loops are located on the corners of a plaquette, with two inequivalent 
configurations 


\ = A (b0) 

a p ~ A [(l,0,0),(0,!)],[(!, 1,0),(1,0)],[(0,1,0),(0,1)] 


^ N t {N t - 1 )u 2Nt+2 


\t 1 , 0 ) 

A [(l, 0,0), (1,0)], [(1,1,0), (0,1)], [(0,1,0), (0,1)] 


= N t u 2Nt {v 4 + w 4 ) 


(4.7) 

(4.8) 


The second contribution is beyond the maximal order considered in this work and neglected. 
It is obvious at this point that there are many different configurations for the four-point 
interaction in addition to the above plaquette form. These appear at u 2Nr+n for n > 4 
and are also neglected here. Once long range interactions of the fundamental and anti¬ 
fundamental loop beyond distance three are included, omission of these terms is no longer 
justified. 

The effective action with the considered additional terms has the following form 

e~ 5 « ff) = e~ s e« (l +\( 2 ^[L 6 ( x )L 6 ( x + rf + L 6 (x)^L 6 (x + i)]J (4.9) 

x,i= 1,...3 

(l + ^(i,i)L8{x)Ls(x + i)^ ( 4 . 10 ) 

x,i= 1,...3 

Y[ (l + X p [L(x)L\x + j)L(x + i + j)L ji (x + j) + c.c.] S j. (4.11) 

x .3, j<i 


We have introduced a short hand notation for the sextet (+6( x ) = X( 2 ,o)(FW( x ))) and 
adjoint (+8( x ) = X(i,i)(kF( x ))) Polyakov loop. 
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4.2 Numerical calculation of the effective couplings 



u(/3) 


Figure 9. The coupling constants of the adjoint next neighbour, sextet next neighbour, and 
plaquette interactions in the effective theory (4.11). These constants are a solution of the equations 
(4.13). The lines are the strong coupling results (4.1), (4.2), and (4.7). The vertical line indicates 
p c of full Yang-Mills theory. The simulations have been done on a 4 x 24 3 lattice. Results for Ai, 
A 2 , and A 5 of the fundamental loop are shown for comparison. 

The additional couplings are related to the correlators of Polyakov loops in higher rep¬ 
resentations and the four-point correlator of fundamental Polyakov loops. The coefficients 
A that follow directly from these measurements are 

A(i,i) = <X(i,i)(V^((0,0,0))) X (i,i)(VlY((l, 0,0)))> 

A (2 ,o) = (X(2,o)(^((0, 0, 0))) X (o,2)(W((l, 0, 0)))) (4.12) 

\ = (L((0,0,0))Lt((l,0,0))L((l,l,0))Lt((0,l,0))) . 

As in the previous case we expand these expectation values in the effective couplings \ n . 
Along the same lines as for the couplings Ai to A 5 additional equations can be derived that 
enlarge the system of equations (3.11). The correlators of the adjoint and sextet repres¬ 
entation receive contributions from the fundamental interactions only at high orders. In 
contrast, the plaquette correlator contains the low order two-point fundamental interaction 
in its disconnected part. At the current order in our investigations we get the following set 
of additional equations 

A(ii) = A( 1;1 ) + 8Af + 16AjA2 + • • • 

A(2 ,o) = A ( 2,0) + 4Af + 8A^A2 + • • ■ (4-13) 

Ap = Xp + 2Aj + 12\f + 24A^A2 + ■ ■ • • 
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The numerical solution for the additional couplings from simulations at N T = 4 is shown 
in Figure 9. In the region of small /3 the couplings of the adjoint and sextet representation 
as well as the four-point interaction are larger than the long distance couplings in the fun¬ 
damental representation. According to our estimates, this is not compensated by the larger 
number of neighbours for the fundamental Polyakov loop interaction at larger distances. 
On the other hand, we observe a stronger increase of the fundamental couplings as (5 gets 
larger and a change of the ordering of the interactions towards the phase transition. 

5 Conclusions 

We have presented a new numerical method to determine the couplings of three-dimensional 
effective Polyakov loop lattice theories for SU(1V) Yang-Mills systems at finite temperature, 
which arise by integrating out the spatial degrees of freedom. The method is based on the 
use of two different character expansions. In the first, the effective couplings are identical 
to definite ?r-point functions of Polyakov loops which can be easily simulated in full Yang- 
Mills theory. While the corresponding truncated effective theory has (within statistical 
accuracy) exact couplings, it will reproduce correlation functions only up to a maximal 
distance, depending on the truncation, beyond which all n-point functions vanish. A 
second character expansion, commonly used in strong coupling approaches, contains long 
range correlations in any truncation, but its couplings are not directly sinrulable. Both sets 
of effective couplings are small in the entire parameter range and related by a perturbative 
expansion with high accuracy, where the ordering of couplings is guided by strong coupling 
power counting. The resulting effective theory has the same structure as the one determined 
in a strong coupling expansion [7], but with non-perturbatively improved couplings which 
can also be determined in the deconfined phase. A large number of effective couplings can 
be extracted by simulations with relative ease compared to an analytic strong coupling 
series, but of course this has to be redone for every set of parameter values (/3,N T ). 

We find the strong coupling series for the nearest neighbour interaction to agree with 
the improved coupling up to the immediate neighbourhood of the phase transition, which 
marks the convergence radius of the strong coupling series. As already observed in [12], 
a quantitative description of the Polyakov loop correlator requires the non-local couplings 
up to the lattice distance of interest. Including these with our new method, quantitative 
agreement is achieved. Also near the phase transition, we observe an increasing relevance 
of long-range couplings, which can be related to the fact that the correlation length is 
maximal there. For N T = 4 and N T = 6 the relevant contributions are the nearest and 
next-to-nearest neighbour interactions, again in accord with earlier observations [12]. As 
the lattice gets finer, more interactions have to be included, depending on the desired 
accuracy. The critical coupling can be obtained at 10% accuracy from the one-coupling 
theory with the analytic effective coupling in a range N r = 2 — 10. Conversely, with our 
numerical determination of the effective couplings, effective theories reproduce the location 
of the phase transition to 1% accuracy with one coupling for N r = 4 and five couplings at 
N r = 8. It would be very interesting to extend this approach to the fermionic sector of 
QCD. 
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